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The Thermal Renormalization Group can be employed to 
study the dynamics of T / Quantum Field Theories close 
to second order phase transitions, where neither resummed 
perturbation theory nor first principle lattice simulations can 
be employed. 

As an example, I discuss the computation of the plasmon 
damping rate in the scalar field theory from T S> Tc down to 
Tc- As the critical point is approached, the lifetime of long 
wavelength thermal fluctuations diverges. 

Taking this effect into account, the notion of quasiparticle, 
and kinetic approaches, can be extended close to the critical 
regime. The consequences on the scenarios for topological 
defect formation in the early universe are also discussed. 



I. CRITICAL SPEEDING UP VS. CRITICAL 
SLOWING DOWN 

The last few years have seen a dramatic progress in 
the understanding of the static properties of relativis- 
tic quantum field theory (QFT) at high temperature. 
Perturbative computations, lattice simulations, and the 
Wilson Renormalization Group, have been successfully 
applied to many different problems, e.g. to the inves- 
tigation of the free energy in the scalar theory, in the 
electroweak standard model, and in QCD. On the other 
hand, there are many interesting dynamical phenomena, 
like the generation of the cosmological baryon asymme- 
try or the formation of topological defects in a second 
order phase transition, in which time plays a crucial role 
and a real time formulation (either in or out of equilib- 
rium) is mandatory. Since first principle lattice simula- 
tions in the real time formalism are not feasible and per- 
turbation theory is plagued by severe infrared (IR) prob- 
lems, semi-phenomenological approaches, mainly based 
on Langevin-type equations, have been employed so far. 

The Thermal Renormalization Group (TRG) intro- 
duced in |l ] offers a promising alternative to this strategy. 
First of all, it is formulated in the real-time formalism, 
where non-static quantities are more conveniently com- 
puted. Second, IR divergences are kept under control by 
the introduction of an explicit momentum cut-off. Third, 
by lowering the cut-off, the IR regime is smoothly ap- 
proached via a continuos sequence of effective theories 
describing the dynamics of the soft modes only, the hard 
ones having been integrated out. This makes the TRG 
particularly suited to the investigation of those dynami- 
cal systems which are dominated by the soft modes. 

In this talk I will discuss the computation of the plas- 
mon damping rate in the real scalar theory Xcj) 4 /4!, from 



temperatures T 3> Tc down to the second order phase 
transition for T — > Tc- The plasmon damping rate is 
defined as 



(rr \ - n/(w fe ,fe) 
lkiT) = 2u k 



(1) 



where oj 2 . = k 2 + m(T) 2 , m(T) is the thermal mass and 
11/ is the imaginary part of the self-energy. 

The physical interpretation of the dynamical quantity 
7k(T) was given by Weldon long ago (see also 0): 
if the plasma is slightly out of thermal equilibrium then 
7fc(T) gives half the relaxation rate of the quasiparticle 
distribution function to its equilibrium value, 



dSrik 
dt 
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(2) 



where Srik is the deviation of the distribution function 
from equilibrium 5n k = Tik — n c k q . 

In the following, I will be interested in the behavior of 
long wavelength fluctuations as the critical temperature 
is approached, so I will concentrate on the damping rate 
for vanishing spatial momentum, "f(T) = 7fe=o(T). 

In perturbation theory, j(T) is a two-loop effect, and 
is given by 
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l/m(T) = 



-1/2 



is the Compton wave- 



where I 
length. 

In refs. || it was shown that the above two-loop result 
can be reproduced in the classical theory provided that 
the Compton wavelength l qu is identified with the classi- 
cal correlation length l c . This can be understood realiz- 
ing that 7 probes the theory at scales u = m(T) <C T (if 
A is perturbatively small), where the Bose-Einstein dis- 
tribution function is approximated by its classical limit, 



N(ui) 



(3u> 



(3m(T) 



(if f3u> « 1) , (4) 



with /3 = 1/T. Considering for instance the 1-1 compo- 
nent of the propagator in the real-time formalism, 
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k 2 
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N(\ko\)), (5) 



we see that when the loop momenta are k <^T the 'sta- 
tistical' contribution to the imaginary part of the propa- 
gator dominates over the 'quantum' one, i.e. N 1/2, 
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and the leading order result can be obtained neglecting 
the T = quantum contributions to the loop corrections. 

The above argument has been employed to motivate 
the use of classical equations of motion in the study of the 
evolution of long wavelength modes in scalar and gauge 
theories. More recently, this approach has been improved 
by many authors including the effect of Hard Thermal 
Loops in the equations of motion for the 'soft' modes (see 
for instance 0j8|). The separation between hard and soft 
modes has been made explicit in ref. || by introducing a 
cutoff A such that m(T) < A < T. 

In the T » T c limit we have m(T)/T ~ (A/24) 1/2 , 
therefore the 'statistical' dominance is valid only up to 
©(A 1 / 2 ) corrections. On the other hand the critical re- 
gion corresponds, by definition, to m(T) — > 0, and the 
statistical limit is exact in this case. 

What does perturbation theory predict in this limit? If 
we trust the two- loop/classical expression in eq. (j^), then 
we see that 7 diverges as the correlation length. Phys- 
ically, this would mean that the lifetime of long wave- 
length fluctuations gets shorter and shorter as the criti- 
cal temperature is approached. However, it is well known 
that (resummed) perturbation theory cannot be trusted 
close to the critical point, due to the divergence of its 
effective expansion parameter, i.e. XT/m(T). 

A more reliable indication of what's going on can be 
obtained from the phcnomenological theory of dynamical 
critical phenomena |) 10 developed by Landau in the 
fifties. Let's consider a non-relativistic, classical scalar 
theory described by an order parameter rj(t, r) and the 
free energy 

T{ji) = J d 3 r[l(dir)) 2 + a(T-Tc)v 2 ] ■ 

If rj is slightly displaced from the minimum of T (at rj = 
0), then it is reasonable to assume that equilibrium will 
be restored at a rate given by 



dy(t,r) 
dt 



= -r 



5T 



Sn(t,r) 



(6) 



where L is a phenomenological parameter which was as- 
sumed to be temperature independent in the original for- 
mulation of the theory, and £(t, r) is a white noise term. 
Eq. @ is usually quoted in the literature as the time- 
dependent Ginzburg-Landau equation. Taking Fourier 
transform, we see that the k = mode vanishes as 
exp(— t/ro) whith the space-independent relaxation rate 
given by 



r = 2 7 (T) - 2aF(T - T c ), 



(7) 



thus vanishing as T — > Tc, in open contrast with the 
divergent perturbative result. 

Dynamical critical phenomena have been the subject of 
intense study by the condensed matter community in the 
seventies (see the reviews in |i7J ) . The Ginzburg-Landau 
approach has been improved by renormalization group 



methods coupled to e or l/N expansions. In general, the 
assumption that transport and kinetic coefficients such 
as F in eq. (^) are temperature independent turns out to 
be false. In any known case, however, the combination 
(T — Tc)T still vanishes. Then, the critical slowing down 
found in (Q) persists, even though with a different power 
law. 



II. WHAT ABOUT 3+1 DIMENSIONAL, 
RELATIVISTIC, TV 0, QFT'S? 

We have seen that the behavior of 7 predicted by per- 
turbation theory is exactly opposite to what is obtained 
by the Landau-Ginzburg approach. On the other hand, 
the latter is based on quite general and physically rea- 
sonable assumptions, while it is well known that pertur- 
bation theory is unreliable in the critical region. 

In refs. [ pTp^ ,[l| it was shown that the key effect which 
is missed by perturbation theory is the dramatic thermal 
renormalization of the coupling constant, which vanishes 
in the critical region. This is just a consequence of trivi- 
ality of A0 4 theory, which implies that, in absence of any 
physical IR cut-off, the renormalized coupling at zero ex- 
ternal momentum vanishes, even if the ultraviolet cut-off, 
Auv, is kept fixed. This is true even at T — 0, where A 
vanishes as 1/ log(m/A^y) when m — + 0. At T ^ 0, due 
to the effective three dimensionality of the theory in the 
critical region, we have a stronger, linear, dependence, 
A(T)~ m (T). 

In the framework of Wilson Renormalization Group 
this can also be understood as follows. The IR regime of 
the four-dimensional field theory at T — Tc is related to 
that of the three-dimensional theory at T = 0. In partic- 
ular, the three-dimensional runningcoupling is obtained 
from the four-dimensional one by |12| 



A 3D (A) = A(A) 



T 

A ' 



where A is the running parameter. At the critical point, 
Asn flows in the IR to the Wilson-Fischer fixed point 
value Ag^, ^ 0, so that the four-dimensional coupling 
vanishes, 



A(A) 



^3D 



(for T ~ T c and A -> 0) . 



The critical exponent governing the vanishing of A is the 
same as that for m(T), so that the ratio X(T)T/m(T) 
goes to a finite value at Tc, and a second order phase 
transition is correctly reproduced. 

We will see that the running of the coupling constant 
for T ~ Tc is crucial also in turning the divergent be- 
havior of eq. (||) into a vanishing one. As a first rough 
ansatz one could just replace the tree coupling A with 
the thermally renormalized one in cq. and readily 
see that the critical speeding up is indeed turned into 
the expected slowing down. However, as we will see, this 
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gives only a qualitatively correct answer, due to the log- 
arithmic singularity of the on-shell imaginary part when 
m(T) vanishes. 

III. THE TRG 

Let me now present in some detail the computation of 
7(T) in the TRG framework. This method allows a com- 
putation of the damping rate for any value of T, from 
very high values, where perturbation theory works, to 
the critical region, where renormalization group methods 
like those of |10[| are necessary to resum infrared diver- 
gencies. Remarkably, the TRG is applicable also in the 
intermediate region, in which none of these methods can 
be employed. 

Wilson's RG idea was originally formulated on the lat- 
tice [ pH , and then implemented in continuum QFT's by 
Polchinski In this framework, it is based on the in- 
troduction of an explicit IR cut-off, via the modification 
of the tree level propagator. Let me briefly recall how it 
works for the scalar theory at T — 0. The starting point 
is the substitution 

p z + 

(where O can be Heavyside's step function or a smooth 
cut-off) in the usual expression for the generating func- 
tional, 

Z A [J}= yp0exp[-(I^.^ 1 .0 + 5M + J.</>)] . 

Z\ [J] generates Green functions in which only the modes 
with p > A have been integrated out. Deriving with 
respect to A the generating flow equation is obtained 

The exact equation (||) provides a non perturbative def- 
inition of the T = QFT. It interpolates between the 
'bare' theory, defined at some ultraviolet scale A = h-uv, 
and the renormalized one at A = 0, in which all quantum 
fluctuations with momenta < p < Ajjy are included. 

A serious problem of this kind of approach is encoun- 
tered when gauge theories are considered. Indeed, the 
momentum cut-off explicitly breaks gauge invariance, 
leading to A-dependent modified Slavnov- Taylor identi- 
ties (ST) which only in the A = limit recover their usual 
form. Even though in principle one can work with this 
new form of the ST's, it is in practice quite non trivial. 
Moreover, no physical meaning can be assigned to the 
theory at non-zero value of A, and the interpretation of 
A as a physical resolution scale is not possible at all. 

When it comes to T ^ QFT's, two possibilities are on 
the ground. One is to work in the imaginary time formal- 
ism and straightforwardly apply the same steps outlined 



above, getting flow equations which interpolate between 
the T — bare theory and the T/0 equilibrium QFT 
]l2| , p~5| | . As long as one is interested in static quantities in 
the critical regime of non-gauge theories, this approach 
works very well. 

The other possibility is the TRG @. The TRG flow has 
three unique characteristics: i) it is formulated in the real 
time; ii) it interpolates between the T = renormalized 
QFT and the T ^ QFT in thermal equilibrium; in) it 
is explicitly gauge invariant JT6|. All this is achieved by 
exploiting the fact that the tree level propagator in the 
real-time formalism (eq. @) is made up by the sum of 
two terms, 

D n (p) = - + 2^(p 2 -m 2 )^(|p |) (9) 

the first one being just the T = propagator and the 
second containing all the statistical information. The 
cut-off is then introduced in the thermal sector only, by 
modifying the Bosc-Einstcin distribution function as 

N(ko) — » NA(ko) = N(k )9(\k\ - A) . (10) 

Now, as A — > oo the modified propagator reduces to that 
of the T = theory, and the generating functional gives 
the Green functions of the T = renormalized QFT 
(i.e. with all the T — quantum fluctuations integrated 
out). By lowering A, the thermal modes of momenta 
\k\ > A are progressively integrated out. Eventually, 
in the A — > limit, the T / theory in equilibrium 
is reached. Moreover, since AT(|po|) comes along with a 
delta function forcing the momentum on the mass shell, 
our cut-off procedure does not break gauge invariance 




IV. COMPUTATION OF 7 (T) 

The flow equations relevant to the computation of the 
plasmon damping rate are schematically reproduced in 
Fig. 1. The RHS's are formally one-loop expression, but 
they contain exact and A-dependent propagators and ver- 
tices, represented by a black dot. The empty dot repre- 
sents the kernel of the evolution equation, which substi- 
tutes the full propagator in the corresponding leg. It is 
given by 

K\,ij(k) = -p A (*MAd)Aa(|fc| -A)JV(|*o|).B0 

where ph{k) is the full, A-dependent, spectral function, 
e(x) — 6(x) — 9(—x), and = 1 with i,j = 1,2 the 
thermal indices. 

We need to compute the real and imaginary parts of 
the (1-1) component of the self-energy 

Sn(w±ie,fc;A) = Ef x (to, k; A) ± iS^ (w, k; A) (11) 

where the real part is the same as that of the self-energy 
appearing in the propagator, Ef 1 (fc;A) = n^(A;; A) and 
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the imaginary part is linked to 11/ in ([!]) by 11/ 



dA 



A. 



dA 



■XX -XX -XX 



FIG. 1. Schematic representation of the evolution equa- 
tions for the two-and four- point functions. 

In order to determine the A-dependent four-point func- 
tion we must consider the system containing also the sec- 
ond equation in Fig. 1 . Since the theory has a - possibly 
spontaneously broken - Zi symmetry, the trilinear cou- 
plings are not arbitrary and will be determined from the 
mass and the quartic coupling as we will indicate below. 
All n-point vertices with n>4 will be neglected. 

The initial conditions for the evolution equations are 
given at a scale A = Ao 3> T (due to the exponential 
damping in the Bose-Einstein function, A 10 T will 
be enough). Here, the effective action of the theory is 
approximated by 



r Ao (<i>) 



where /i\ o and Aa are the renormalized parameters of 
the T = theory. We are interested to the case in which 
the ^-symmetry is broken at T = 0, so we will take 
Ma < 0. 

Then, we make the following approximations to the 
full propagator and vertices appearing on the RHS of the 
evolution equations: 

i) the self-energy in the propagator is approximated by a 



running mass, II(fc; A) 



given by 



ml = i? K (if a4 > 0) , -2/x A (if fil < 0) 
ii) the four-point function is approximated as 

Hi) the three-point function is given by 



r 



(3) 



(p,) ~ (if a4 > 0) , ^-GAamI (if A < 0) 



Moreover, all vertices with at least one of the thermal 
indices different from 1 will be neglected. 

Some comments are in order at this point. In pertur- 
bation theory, the on-shell imaginary part is given by the 
two-loop setting sun diagram. The evolution equations 
of Fig.l contain instead only one-loop integrals, so how 
can an imaginary part emerge? The crucial point here 
is that the momentum dependence of the imaginary part 
of the four-point function has to be taken into account 



(see point ii)). When the latter is inserted into the upper 
equation, 11/ is generated. 

As long as the A-dependent action is in the broken 
phase, trilinear couplings will be present. By neglecting 
the imaginary part of the self-energy in the propagators 
on the RHS (point i) above), we loose the contribution to 
Hi obtained by cutting the second diagram in the RHS 
of the upper equation of Fig.l. However, if we restrict 
ourselves to temperatures T > Tq, the trilinear couplings 
will be different form zero only in a limited range of A. 
They will not contribute in the IR, where the dominant 
contributions to the on-shell imaginary part emerge (see 
Fig. 3). The error induced by this approximation can be 
estimated noticing that the contribution we are neglect- 
ing is of the same nature as the one obtained by inserting 
the last diagram of the second equation in the equation 
for the two-point function. The latter is taken into ac- 
count and its effect on the full imaginary part is of the 
order of a few percent. 

We have now a system of four evolution equations for 



n*(A;A) 



'A- 



Sf 1 (A;;A), A A , and r] A {Pi), with ini- 



tial conditions m Ao = — 2/i Ao , Aa q , and S{ 1 (/c;Ao) 
??Ao (Pi) = 0; respectively. We give, for simplicity, the ex- 
plicit form of the evolution equations in the symmetric 
phase only, where the trilinear coupling vanishes; 

8 a A 3 N(w A ) 

A— m A = ~—r A A 

dA 47r 2 iv A 



8 A 3 d f N(u A ) \ a 

A dA AA - ^I^dmJ J Aa 

OA 47T Z U1 A 

d A 2 
A-^q A {k, -k, q A , -q A ) = — A(C A (fc Q - uj Aj k) + 

C A (k + cu A ,k)) 

where uj a = (A 2 + m^) 1 / 2 , k = (fc ,0), q A = (cj a , \q\ 
A), and 



C A (q ,q) = L^IA. (J (go ) + J(_ 9o) ) ; 
87T lo a \q\ 



with 



J (go) 



if 



2 8((A 2 + ql - 2u; A qo) 1/2 
2uj A q a -ql + \q\ 2 



A) 



2|q|A 



< 1 



otherwise. 



Notice that the subsystem for m 2 and Aa is closed and 
can be integrated separately. 

We then proceed as follows. Fixing the temperature, 
we first integrate the subsystem for m 2 and Aa down to 
A = in order to find the plasmon mass, Wa=o- Then 
we fix the external momentum of Hi on this mass-shell, 
k — (m A=0 , 0) and integrate the full system from A = Ao 
down to A = 0. 
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t=(T-T c )/T c 



FIG. 2. Temperature dependence of the coupling constant 
(solid line), and of the damping rate with the effect of the 
running of A included (dash-dotted) and excluded (dashed). 
The values for 7 have been multiplied by a factor of 10. 

In Fig. 2 we plot the results for the damping rate 7 
and the coup ling constant at A = 0, as a function of the 
temperature [ p.8| . The dashed line has been obtained by 
keeping the coupling constant fixed (A-indepcndent) to 
its T — value (A = 10~ 2 ), and reproduces the divergent 
behavior found in perturbation theory (eq. (^)). The cru- 
cial effect of the running of the coupling constant is seen 
in the behavior of the dot-dashed line. For temperatures 
close enough to Tc, the coupling constant (solid line in 
Fig. 2) is dramatically renormalized and it decreases as 

A A=0 (T) ~ f 

where t = (T - T c )/T c and we find v ~ 0.53 Q. The 
mass also vanishes with the same critical index. The de- 
creasing of A drives 7 to zero, but with a different scaling 
law, 

7A=oCr)~f / io g t. 



1 More accurate results for the critical exponents from the 
TRG have been obtained recently in ref. []19|. 



The above expression can be understood noticing that 
the two-loop contribution to II/, computed at vanishing 
u; = m(T), goes as A 2 log m(T). The RG result replaces 
A with the renormalized coupling, then from (Q) we have 

A 2 _ 

7~ -^=2-logm A ^o~*"log* ■ 
to a=0 

Taking couplings bigger than the one used in this let- 
ter (A = 10~ 2 ), the deviation from the perturbative 
regime starts to be effective farther from Tc- Defin- 
ing an effective temperature as Aa = o(T)/Aa < 1/2 for 
T c < T < T e ff we find that t e f-f scales roughly as 
teff ~ A A[) . 

In Fig. 3 we plot the running of Aa(T) and Ja{T) for 
two different values of the temperature. When T ^> Tc 
most of the running takes place for A <; A so ft, so it is 
safe to stop the running at this scale neglecting the effect 
of soft loops. When T — > Tc this is not possible any 
longer, since most of the running of A and 7 takes place 
for A < A so f t . Thus, the vanishing of the mass gap forces 
us to take soft thermal momenta into account. 



X(t=20) 




11N. \ " 

w' 7 w' 6 10'* io~ l 10' 3 w' 2 10'' 1 10 

A/T 



FIG. 3. A-running of the coupling constant A and of the 
damping rate 7 for two different values of t = (T — Tc)/Tc- 
The values of 7 (dot-dashed lines) have been multiplied by 
500. 



V. CONCLUSION AND CONSEQUENCES 

The world of dynamical critical phenomena in rela- 
tivistic QFT's at T ^ is still largely unexplored, the 
reason being that the two main tools employed in the 
static case, i.e resummed perturbation theory and lat- 
tice simulation, are of no use in this context. I think 
that the TRG is on a much better shape, and can be 
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properly employed to study time dependent correlation 
functions, relaxation rates, and transport coefficients at 
any temperature, from the critical one to much higher 
(or lower) ones. 

In the case of the plasmon damping rate discussed here, 
the expected critical slowing down of long wavelength 
fluctuations has been obtained from first principles. The 
dominant effect has been shown to be the thermal renor- 
malization of the coupling constant, which turns the di- 
vergent behavior of the plasmon damping rate into a van- 
ishing one. In the future this work will be extended to 
non-zero spatial momentum, and to higher orders of ap- 
proximation of the flow equations in order to check the 
scaling relations between dynamical critical exponents in 
a non-trivial way |ic|] . 

Before concluding, let me remark two interesting phys- 
ical consequences of this result. 



FT 



'Yno run/ 



r y/m 



possible to talk about quasiparticles around the critical 
point. In particular, no kinetic description close to the 
phase transition could be conceivable. The TRG result 
shows that it is not necessarily the case, as we see in 
Fig. 4. Now 7/m diverges only as 1/logi, and it may 
be less then unity very much closer to the critical point. 
Indeed, for A = 10~ 2 , the ratio j/m(T) becomes larger 
than unity for t <, 10 -5 in perturbation theory, whereas 
the RG result is still j/m(T) ~ 0.3 at t ~ 10" 9 . 

The second consequence concerns second order phase 
transitions in the early universe. In a cosmological set- 
ting, the increasing lifetime of the fluctuations of the or- 
der parameter may modify the dynamics of second order 
- or weakly first order - phase transitions. Indeed, as the 
critical temperature is approached, the expansion rate 
of the universe exceeds the thermalization rate of long 
wavelength fluctuations, and the latter freeze out. There 
will be a range of temperatures around Tc in which the 
equilibrium expressions for the effective potential or the 
fluctuation rates cannot be employed. In brief, second 
order phase transitions in the early universe take place 
out of equilibrium. 

At the electroweak scale this has no practical conse- 
quences. Due to the smallness of the expansion rate at 
that epoch, the temperature interval in which there is de- 
parture from equilibrium is too short to have any effect. 

On the other hand, at higher temperatures, typically 
for T J; 10 10 GeV, the effect may be very strong and must 
be taken into account. This would lead to a scenario for 
the formation of topological defects similar to that pro- 
posed by Zurek in |p0|, and discussed by J. Rivers at this 
conference. In this scenario, the initial length scale of 
the network of topological defects is given by the correla- 
tion length at the freeze-out, £/. ., whose determination 
requires in turn the knowledge of the temperature depen- 
dence of j(T), computed here. 
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t=(T-T c )/T c 



FIG. 4. The ratio 7/m vs. T as obtained from TRG 
(dot-dashed line) and without including the running of A 
(solid line) 

The first consequence has to do with the notion of 
quasiparticle. In order for this concept to make sense, 
a narrow resonance is needed, a good criteria being 
7(T)/m(T) Ss 1. In perturbation theory this quantity 
diverges as t~ v , and one would conclude that it is not 
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